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The analyzing power A y for the pp —> ppr\ reaction has been determined at the beam mo- 
mentum pbeam = 2010 MeV/c, corresponding to the excess energy Q = 10 MeV. In the 
paper the method of the data analysis is briefly presented. 
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1 Introduction 

The close-to-threshold measurements of the total cross section for the pp — * ppr\ reaction [1-6], 
the measurements of the differential cross sections [7-9], the high statistic investigations of the 
pp?7 final state dynamics [ 1 0, 1 1 ] as well as the results of the first ever analyzing power experiment 
for the pp — ► ppr\ reaction [12] are in agreement with the hypothesis that the main contribution to 
the production of the r\ meson in the proton-proton collisions comes from the mesonic excitation 
of the Sn(1535) resonance and it's further decay into the p?; pair [13, 14], as depicted in Fig.^ 
However, it is still not settled what are the relative contributions to the production mechanism 
originating from each meson exchange denoted in Fig. ^ A possible way for elucidation of 
this problem is a precise measurement of the polarization observables, as they are very sensitive 
to the model assumptions on the type of the meson being exchanged in order to excite the Sn 
resonance. For example, from Fig. where the predictions of the analyzing power for the 
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pp — > pp?7 reaction at the excess energy Q = 10 MeV are shown, it can be noticed that the 
analyzing power function strongly depends on whether we deal with pseudoscalar [15] (full line) 
or vector meson exchange model [13] (dotted line). 
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Fig. 1. The Feynman diagram for one of the possi- 
ble explanations of the rj meson production in ele- 
mentary proton-proton collisions. ISI and FSI stand 
for the initial and final state interactions, respec- 
tively. 



Fig. 2. Predictions of the analyzing power for the 
pp —* ppr\ reaction at Q=10 MeV as a function of 
the center-of-mass (CM) polar angle of r\ (6cm)- 
Full line are the predictions of the pseudoscalar me- 
son exchange model [15], whereas the dotted line 
represents the results of the calculations, based on 
the vector meson exchange [13]. 



Therefore, in order to verify the validity of these two models, an experiment devoted to a 
determination of the analyzing power for the pp — » ppr\ reaction has been performed at the stor- 
age ring COSY in the Research Center Jiilich in Germany. In this article we present preliminary 
results obtained at excess energy Q = 10 MeV, where the predictions of the models differ at most. 



2 Experimental method 

In the measurements a polarized proton beam [16,17] with the momentum of Pbeam = 2010MeV/c 
has been used. This beam momentum value corresponds to the excess energy Q = 10 MeV. The 
beam of protons has been scattered on the H2 clusters, which were produced by the hydrogen 
cluster target [18]. The COSY-11 detection setup [19,20] has been used for the detection of the 
reaction products. 

As the COSY- 1 1 system is a one-arm detection setup, in order to register both types of events: 
scattered to the right and to the left side (with respect to the polarization plane) there is a need 
to perform the experiment with the separate cycles, in which the spin of the proton beam has 
to be flipped, as schematically shown in Fig. [5] In this figure we also define what is meant by 
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the terms: "scattering to the right" and "to the left". This definition strictly follows the Madison 
convention [21]. 



Pin 



a) 



R(sp±n "up") 



b) 



L (spin "down ") 



Fig. 3. Definition of the scattering to the right (a) and to the left side (b) with respect to the polarization 
plane. P is a polarization vector of the beam of protons, whereas pi n and p out denote the momenta of 
incoming and outgoing particles, respectively. The polarization vector of Fig.|3|a points out of the paper, 
whereas the polarization vector of Fig.|5|b points into the paper. 



In the following we define the polarization plane as a plane spanned by the beam momentum 
vector pbeam = [0, 0, pl eam \ and a polarization vector P = [0, P, 0]. If the incident beam 
consists of spin ~ particles with a transverse polarization, which was the case in our experiment, 
the formula for the differential cross section a(9, <fr) for a scattering into the solid angle around 
9 and <\> angles reads: 



a (9)(l + A v (9)P-n), 



(1) 



where <7o (9) denotes the differential cross section for a scattering of an unpolarized beam, A y (9) 
is an analyzing power of the reaction, and n is a unit vector along p. m x p ou t- <j> is an angle 
between P and h, i.e. P ■n = P cos <j> sin 9 and 9 denotes the angle between pi n and p ou t- In the 
experiment we were restricted to the plane corresponding to <p = 0° in case of spin "down" (or 
4> = 180° for spin "up"), where the detection efficiency of the COSY-1 1 setup is the highest and 
it decreases drastically when going outside this particular scattering plane (see Fig.|4jl. Therefore, 
in what follows we will consider the scattering in this plane solely. 
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Fig. 4. Monte-Carlo simulation of the distribution of the events as a function of the angle . 



The total production rate to the right N R and to the left side N l can be expressed as follows: 

N H (0)=N(9,4> = K)=(T R (0,ir)E R {0,n) J L R dt R = a (6)(l-A y (6)P R )E R (6,ir) J L R dt R , 

(2) 



N L (6) = N(6,cj> = 0)=a L (e,0)E L (6,0) J L L dt L = <T O (0)(l+A y (6)P L )E L (6, 0) J L L dt L , 

(3) 

where the (9, tt) = El (0, 0) is the function of the efficiency of the COSY- 1 1 detection system, 
L R / L stands for the luminosity for spin up/down, and t R /i, denotes the time of the measurement 
with spin up/down, respectively. The polarization degrees for spin up and for spin down 
are equal within the 5% of accuracy, as has been shown in the previous measurements by means 
of the EDDA polarimeter [22]. Therefore for the further considerations we will assume that 
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Pl=Pr = P. Dividing equation[3]by|2l and introducing a relative luminosity: 

J Ll^l 

results in the following formula for the analyzing power: 

A f cos 9)-- L reiN L (cose)-N R (cos6) 
y[ j " PL rel N L ( C os9)+N R (co S ey W 

Therefore, the values of the relative luminosity L re ;, the spin-averaged polarization degree 
P, together with the production rates to the right Nr, and to the left side are the physical 
observables to be found in order to determine the analyzing power function A y . In the following 
subsections the methods used to obtain these quantities will be presented. 



2.1 Calculations of the relative luminosity and polarization 

The experiment has been performed with the 300 s long cycles, for which the spin of the incident 
proton beam has been flipped from cycle to cycle. This method was intended to reduce the 
systematical uncertainties caused by the change of the beam parameters due to the variation of 
the target densities. The duration of the cycle has been set up in such a manner, that it was 
significantly shorter than the time scale (10 hours) for the substantial changes of the density of 
the target. 

For the determination of the relative luminosity a detection system consisting of two pairs of 
the scintillators placed in the polarization plane have been used [23] . Due to the parity invariance, 
a cross section for any scattering into this plane does not depend on the polarization degree. Thus, 
the number of the coincidences n between scintillators has been used as a measure of the absolute 
luminosity. Following the definition given in eq.|4] we can express the relative luminosity as a 
ratio of the number of coincidences in the polarization plane during the cycles with spin "up" 
(tir) and "down" (n/,): 



The relative luminosity determined using eq.|6]was found to be L re ;=0.96468± 0.00065. 

The detailed method of determination of the averaged polarization value has been described 
in [23]. For the calculations the assymetry of the elastically scattered events in the plane per- 
pendicular to the polarization plane has been found for different scattering angles in the center- 
of-mass system (CM). In order to determine the values of the analyzing powers for the beam 
momentum of Pbeam = 2010 MeV, the linear interpolations between available experimental re- 
sults for the proton-proton elastic scattering at Pbeam — 1995 and 2025 MeV/c have been used. 
Two data sets of the analyzing powers for the pp — ► pp reaction for the mentioned above beam 
momenta, as measured by the EDDA collaboration [22] at different CM scattering angles, are 
presented in Fig. EI 
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Fig. 5. Analyzing powers for the pp — > pp elastic 
scattering as determined by the EDDA collabora- 
tion [22] at the beam momenta specified inside a 
figure. To obtain the analyzing powers for a beam 
momentum used in the experiment a linear approx- 
imation has been applied. 



Fig. 6. Average polarization degree versus the time 
of the measurement. Open circles are the results ob- 
tained by means of the COSY polarimeter, whereas 
the full dots are the polarization values determined 
using the method described in [23]. Data obtained 
using the independent methods are in line. 



Fig-IUdepicts the variation of the beam polarization during the time of the measurement. Data 
sets obtained using the COSY polarimeter [24](open circles) are confronted with the results of 
the measurements with the COSY-1 1 polarimeter [23,25,26](fullcircles). An agreement between 
both sets of data is visible in Fig.|6] The average value of the beam polarization for the whole 
period of the measurement was found to be P=0.680± 0.010. 

2.2 Determination of the background free production rates for the pp — » ppr\ reaction 

A typical missing mass spectrum for the pp — » ppr\ reaction as measured using the COSY-1 1 
detection setup is presented in Fig. [7] Figure shows the spin-averaged missing mass spectrum, 
gathered during the whole time of the experiment. Over the wide multi-pionic background a 
clear peak originating from the r\ meson production is visible. 

Optimizing the statistics and the expected shape of the analyzing power function, the range 
of the center-of-mass polar angles of the r\ meson emission has been divided into four bins. To 
separate the actual production rates from the background, the reactions with multi-pionic produc- 
tion as well as the events with r] production have been simulated, using a program based on the 
GEANT3 [27] code. Generated events which fulfilled conditions equivalent to the experimental 
trigger have been analyzed in the same way as the experimental data. In such a way we have 
obtained the missing mass shapes of the background reactions (pp — > pp2ir, Sir, A-k) as well as 
the shape of the signal (pp — > pprf). In order to perform credible simulations of the missing mass 
spectra the geometry of the drift chambers (used for the momentum reconstruction) as well as the 
position and geometrical parameters of the target have to be known precisely. The angle and the 
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Fig. 7. Spin-averaged missing mass spectrum for the pp 
as measured by means of the COSY- 1 1 detector setup. 



pprj reaction at the excess energy Q = 10 MeV, 



relative position of the COSY- 1 1 drift chambers (2 parameters), a target position (3 parameters) 
and a relative shift between the beam and a target ( 1 parameter) have been optimized using the 
elastically scattered events and the method described in [1 1, 29]. 



The multi-dimensional fit of the simulated missing mass spectra to the corresponding ex- 
perimental histograms has been performed using the MINUIT [28] minimization package. The 
fit has been performed simultaneously to 8 histograms (see Fig. [8}, each with 60 points. The 
X 2 of the fit have been minimized as a function of six parameters: scaling factors of the gener- 
ated background and signal reactions (altogether 4 parameters) and 2 parameters responsible for 
the spread and the absolute value of the beam momentum. The detailed process of the search 
for the best values of these parameters will be described elsewhere [30], here we would only 
like to report that the minimum value of \ 2 was found to be 1 .62, corresponding to the spread 
of the beam momentum equal to 0.2 MeV/c and the shift from the nominal beam momentum 
equal to 2.1 MeV/c. These values are in a good agreement with results of previous COSY- 11 
experiments [29]. 



InFig.|8]the missing masses for the individual 6cm subranges, separately for spin up (upper 
panel) and down (lower panel), are shown. Full dots correspond to the experimental data, the 
dotted line depicts the shape of the background, whereas the solid line represents the best fit of 
the Monte-Carlo data to the experimental spectra. 
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Fig. 8. Missing mass spectra for different ranges of cos 9cm for spin up (two upper panels) and spin down 
(two lower panels). Dots represent the experimental points. Dotted line shows the generated multi-pionic 
background. The full line is a best fit of the sum of the signal and background to the experimental data. 
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3 Results 

Preliminary results of the analyzing power function Aj,(cos 9) are presented in Fig.[9]as the full 
dots. The dotted line shows the predictions determined according to the vector meson exchange 
model [13], whereas the solid line refers to the pseudoscalar meson exchange model [15]. The 
errors are of the statistical nature only. 
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Fig. 9. Preliminary analyzing power function for the pp — > ppr\ reaction at Q=10 MeV. Vertical bars denote 
the statistical error, whereas the horizontal bars stand for the ranges of averaging. 



4 Conclusions 

The preliminary analyzing power results are bared with rather large statistical uncertainties, 
therefore a statement on the mechanism of the 77 meson production relying on this set of data 
is currently unclear. As can be seen from Fig.[9]the data point at cos 6 = 0.75 tend to prefer the 
pseudoscalar meson exchange model [15], whereas the data point at cos# = —0.75 is slightly 
more in line with the predictions of the vector meson exchange model. The results show rather 
small values of the analyzing power in the close-to-threshold region, which may be the indication 
of the 77 production to the s-wave final state, solely. However, for this statement the exact partial 
wave analysis remains to be done. 
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